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Abstract. Infrared (IR) luminosity is fundamental to understanding the cosmic star 
formation history and AGN evolution, since their most intense stages are often obscured 
by dust. Japanese infrared satellite, AKARI, provided unique data sets to probe these 
both at low and high redshifts. The AKARI performed an all sky survey in 6 IR bands 
(9, 18, 65, 90, 140, and l6Qpm) with 3-10 times better sensitivity than IRAS, covering 
the crucial far-IR wavelengths across the peak of the dust emission. Combined with 
a better spatial resolution, AKARI can measure the total infrared luminosity (Ltir) of 
individual galaxies much more precisely, and thus, the total infrared luminosity density 
of the local Universe. In the AKARI NEP deep field, we construct restframe ^pm, 
\2pm, and total infrared (TIR) luminosity functions (LFs) at 0.15< z <2.2 using 4128 
infrared sources. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 
1,9, 11, 15, 18, and 2Apm) by the AKARI satellite allows us to estimate restframe 2>pm 
and \2pm luminosities without using a large extrapolation based on a SED fit, which 
was the largest uncertainty in previous work. By combining these two results, we reveal 
dust-hidden cosmic star formation history and AGN evolution from z=0 to z-2.2, all 
probed by the AKARI satellite. The next generation space infrared telescope, SPICA, 
will revolutionize our view of the infrared Universe with superb sensitivity of the cooled 
3m space telescope. We conclude with our survey proposal and future prospects with 
SPICA. 


1. Lessons from AKARI 
1.1. Background 

Revealing the cosmic star formation history is one of the major goals of observational 
astronomy. However, UV/optical estimation only provides us with a lower limit of 
the star formation rate (SFR) due to obscuration by dust. A straightforward way to 
overcome this problem is to observe in the infrared, which can capture star formation 
activity invisible in the UV. The superb sensitivities of Spitzer and AKARI satellites 
have revolutionized the field. 

Univer se, often used IR LFs are from the IRAS (e.g., ISanders et al. 
201 Ibh from 1980s, with only several hundred galaxies. In addition, 
bolometric infrared luminosities (L/^ 8_iooo//m) of local galaxies were estimated using 
equation in Perault (1987), which was a simple polynomial, obtained assuming a simple 
blackbody and dust emissivity. Furthermore, the reddest filter of IRAS was 100/rm, 
which did not span the peak of the dust emission for most galaxies, leaving a great deal 
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Figure 1. (left) Emission line ratios used to select AGNs from the AKARI all 
sky sample. The contour shows distribution of all galaxies in the SDSS with r < 
\1.11 (regardless of IR d etection). The dotted line is the criterion between starbursts 
and AGNs described in (iKewlev et al.ll200l1) . The dashed line is the criterion by 
(iKauffmann et alJl20Q3t) . Galaxies with line ratios higher than the dotted line are 
regarded as AGNs. Galaxies below the dashed line are regarded as star-forming. 
Galaxies between the dashed and dotted lines are regarded as composites. The blue 
and red dots are for ULIRGs, and LIRGs, respectively. The green squares are Seyfert 
1 galaxies identified by visual inspection of optical spectra. 

Figure 2. (right) Fracti ons of AGN and co mposite galaxies as a function of Lir. 
AGN are classified using (iKewlev et al.ll200ll) among galaxies wit h all 4 lines mea¬ 
sured. Composite galaxies include those classified as AGN using IKauffmann et al.l 

(I2n03h . 


of uncertainty. Using deeper AKARI all sky survey data that cover up to 160yum, we 
aim to measure local L/r, and thereby the IR LF more accurately. 

At higher redshifts, most of the Spitzer work relied on a large extrapolation from 
24/tm flux to estimate the 8, \2pm or total infrared (TIR) luminosity, due to the limited 
number of mid-IR filters. AKARI has continuous filter coverage across the mid-IR 
wavelengths, thus, allowing us to estimate mid-IR luminosity without using a large k- 
correction based on the SED models, eliminating the largest uncertainty in previous 
work. By taking advantage of this, we present the restframe 8, \2pm, and TIR LFs, 
and thereby the cosmic star formation history derived from these using the AKARI 
NEP-Deep data. 

1.2. AKARI All Sky Survey: low-z Universe 

AKARI performed an all-sky survey in two mid-infrared bands (centered on 9 and 18 
pm) and in four far-infrared bands (65, 90, 140, and 160//m). In addition to the much 
improved sensitivity and spatial resolution over its precursor (the IRAS all-sky survey), 
the presence of the 140 and 160//m bands is crucial to measure the peak of the dust 
emission in the EIR wavelength, and thus the L/r of galaxies. We have cross-correlated 
the AKARI EIS bright source catalog with the SDSS DRV galaxy catalog, obtaining 
2357 cross-matched spectroscopic redshifts. 

It is fundamental to separate IR contribution from two different physical processes; 
the star-formation and AGN activity. In Eig. [T] we use [A//]///a vs [OIII]IHp line 
ratios to classify galaxies into AGN or SEG (star-forming galaxies). It is interesting 
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Figure 3. Infrared luminosity function of AKARI-SDSS galaxies. The L/r is 
measured using the AKARI 9, 18, 65, 90, 140 and 160/im fluxes through an SED fit. 
Errors are computed using 150 Monte Carlo simulations, added to a Poisson error. 
The dotted lines sho w the best-ht doub le-power law. The green dotted lines show IR 
LE at z=0.0082 by (iGoto et al.ll201 lb|) . The dashed-dotted lines are higher redshift 
results from the AKARI NEP deep held (iGoto et al.ll2010alf3i . 


that the majority of (U)LIRGs are aligned along the AGN branch of the diagram, im¬ 
plying the AGN fraction is high among (U)LIRGs. This is more clearly seen in Fig. 
m where we plot fractions of A GN as a fun ction of Ljr. These results agree with pre¬ 
vious AGN fraction estimates dGotol l2005h . Improvement in this work is that due to 
much larger statistics, we were able to show fractions of AGN in much finer luminosity 
bins, more accurately quantifying the increase. Especially, a sudden increase of /agn 
at logL/R >11.3 is notable due to the improved statistics in this work. 

For these galaxies, we estimated total IR luminosities (LjA by fitting the AKARI 


photometry with SED templates. We used the LePhare codeU to fit the infrared part 
(>7;um) of the SED. We fit our AKARI FIR photometry with the SED templates from 
Chary & Elbaz (2001; GHEE hereafter), which showed most promising results among 
SED models tested by (iGoto et al.ll2011bh . 

With accurately measured L/r, we are ready to construct IR EEs. Since our sample 
is flux-limited at r = 17.7 and Sgo^,„ = OJJy, we need to correct for a volume effect to 
compute EEs. We used the l/Vmax method. We estimated errors on the EEs with 150 
Monte Carlo simulations, added to a Poisson error. 

In Fig. [5J we show infrared EE of the AKARI-SDSS galaxies. The median redshift 
of our galaxies is z^O.031. 

Once we measured the EE, we can estimate the total infrared luminosity density 
by integrating the EE, weighted by the luminosity. We used the best-fit double-power 
law to integrate outside the luminosity range in which we have data, to obtain estimates 
of the total infrared luminosity density, D/r. Note that outside of the luminosity range 
we have data (Ljr > or Ljr < lO^-^L©), the EEs are merely an extrapolation 

and thus uncertain. 

The resulting total luminosity density is fliR= x 10^ LoMpc“^. Errors are 

estimated by varying the fit within Icr of uncertainty in EEs. Out of O/r, 1.1±0.1% is 


'http://www.cfht.hawaii.eciu/~amouts/lephare.html 
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Figure 4. (left) Restframe Spm LFs. The blue diamonds, purple triangles, red 
squares, and orange crosses show the 8//m LFs at 0.38 < z < 0.58,0.65 < z < 
0.90,1.1 < z < 1.4, and 1.8 < z < 2.2, respectively. The dotted lines show analytical 
fits with a double-power law. Vertical arrows show the 8jum luminosity correspond¬ 
ing to the flux limit at the central redshift in each redshift bin. (middle) Restframe 
12jum LFs. The blue diamonds, purple triangles, and red squares show the 12jum 
LFs at 0.15 < z < 0.35,0.38 < z < 0.62, and 0.84 < z < 1.16, respectively, (right) 
TIR LFs. The redshift bins used are 0.2< z <0.5, 0.5< z <0.8, 0.8< z <1.2, and 
1.2<z <1.6. 



AGN 



Figure 5. (left) Evolution of TIR luminosity density based on TIR LFs (red cir¬ 
cles), 8/rm LFs (stars), and 12pm LFs (filled triangles). The blue open squares and 
orange filled squares are for LIRG and ULIRGs only, also based on our Ltir LFs. 

Figure 6. (right) Evolution of TIR luminosity density by AGN. Results from the 
AKARI all sky survey is plotted with stars at z=0.0082. The red, blue and orange 
points show IR luminosity density from all AGN, from LIRG AGN only, and from 
ULIRG AGN only . Higher redshift results are from the AKARI NEP deep field 
(iGoto et al.ll2010bl) . with contribution from star forming galaxies r emoved. Brown 
triangles are computed from the 8yum LEs (iGoto et al.ll2010bh . 
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produced by LIRG (L/r > lO^^Lf:)), and only 0.03±0.01 % is by ULIRG (L jr > lO'^L©). 
Although these fractions are larger than at z=0.0081 (iGoto et alJl2011b[) . still a small 
fraction of Q/r is produced by luminous infrared galaxies at z=0.031, in contrast to the 
high-redshift Universe. 

1.3. AKARINEP Deep Field: high-z Universe 

The AKARI has observed the NEP deep field (0.4 deg^) in 9 filters {N2, N3, N4, S1,S 9W, 
511.L15.L18 W and L24) to th e depths of 14.2, 11.0, 8.0, 48, 58, 71, 117, 121 and 
275 juJy (5cr, IWada et ^ l2008h . This region is also observed in BVRi'z' (Subaru), 
u' (CFHT), FUV,NUV (GALEX), and J,Ks (KPN02m), with which we computed 
photo-z with ^-0.043. Objects which are better fit with a QSO template are removed 
from the analysis. We used a total of 4128 IR sources down to 100 gJy in the L18 filter. 
We compute EEs using the 1/Emax method. Data are used to 5cr with completeness 
correction. Errors of the EEs are from 1000 realization of Monte Carlo simulation. 


1.3.1. Restframe 8yum luminosity functions 

Monoc hromatic 8/rm lum inosity (L^um) is known to correlate well with the TIR lumi- 


ivtonoc nromatic oum tum mosity is Known to corretate wett witn me i IK tumi- 

nosity dGoto et al.ll2011al) i. especially for star-forming galaxies, because the rest-frame 
8/rm flux is dominated by prominent PAH features such as at 6.2, 7.7 and 8.6 gm. The 
left panel of Fig. |4] shows a strong evolution of 8/rm EEs. Overplotted previous work 
had to rely on SED models to estimate Lg^„, from the Spitzer S 2 Aiim in the MIR wave¬ 
lengths where SED modeling is difficult due to the complicated PAH emissions. Here, 
AKARFs mid-IR bands are advantageous in directly observing redshifted restframe 
8//m flux in one of the AKARFs filters, leading to more reliable measurement of 8-/rm 
EEs without uncertainty from the SED modeling. 


1.3.2. Restframe 12jum luminosity functions 

12/rm luminosity (L n,,m) represents mid-IR conti nuum, and known to correlate closely 
with TIR luminosity dPerez-Gonzalez et al.ll200^ . The middle panel of Fig. |4] shows a 
strong evoltuion of 12yum EEs. Here the agreement with previous work is better because 
(i) 12yum continuum is easier to be modeled, and (ii) the Spitzer also captures restframe 
12/rm in S 24 ^m at z= 1. 


1.3.3. Total infrared luminosity functions 

Easily , we show the TIR EEs in the right panel of Fig. |4l We used lEasache et al. 


d2003l) ’s SED templates to fit the photometry using the AKARI bands at >6/rm {S 7, S 9W, 
S11,L15,L18W and L24). The TIR EEs show a strong evolution compared to local 
EEs. At 0.25 < z < 1.3, Ll,„ evolves as oc (1 -i- 


1.3.4. Cosmic star formation history 

We fit EEs in Fig. |4] with a double-power law, then integrate to estimate total infrared 
lu minosity density at various z- The restframe 8 and 1 2;um EEs are converted to Ltjr us¬ 
ing iPerez-Gonzalez et ^ d2005h : ICaputi et all d2007h before integration. The resulting 
evolution of the TIR density is shown in Fig. |5] The right axis shows the star formation 
density assuming Kennicutt (1998) . We obtain oc (1 -i- z)^'^“°'^. Comparison to 

Q.UV bv ISchiminovich et al.l ( 2005 ) suggests that Qtir explains 70% of fltotai at z=0.25. 
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Figure 7. Comparison of survey speed for 5x5 arcmin^ area (1 hour, S/N=5, 
R=5). SPICA’s confusion limit is shown in the blue filled-circles. 


and that by z^l.3, 90% of the cosmic SFD is explained by the infrared. This implies 
that Q.TIR provides good approximation of the Q-totai at z > 1. 

In Fig. m we also show the contributions to Q-tir from LIRGs and ULIRGs. From 
z=0.35 to z^l.4, Q.JR by LIRGs increases by a factor of ~1.6, and Q//; by ULIRGs 
increases by a factor of ~10. More details are in Goto et al. (2010a). 

1.3.5. Cosmic AGN accretion history 

We have separated the from Therefore, we can also investigate 

By integrating IR LFagw> we show the evolution of in Fig. which shows 

a strong evolution with increasing redshift. At a first glance, both and 

show rapid evolution, suggesting that the correlation between star formation and black 
hole accretion rate continues to hold at higher redshifts, i.e., galaxies and black holes 
seem to be evolving hand in hand. When we fit the evolution with (1+z)^, we find 
oc(l+z)‘^-^-*’'^. A caveat, however, is that estimated in this work is likely to 
include IR emission from host galaxies of AGN, although in optical the AGN compo¬ 
nent dominates. Therefore, the final conclusion must be drawn from a multi-component 
fit based on better sampling in FIR by Herschel or SPICA, to separate AGN/SFG con¬ 
tribution to Ljr. The contribution by ULIRGs quickly increases toward higher redshift; 
By z=L5, it exceeds that from LIRGs. Indeed, we found Q.f^^{ULIRG) oc(l-i-z)^-^-°-® 
and af^^(LIRG) oc(l-rz)5-4±0-5. 


2. Wide area survey with SPICA’s Mid-infrared Camera and Spectrograph (MCS) 

SPICA dNakagawa et al.ll201 ih is the next-generation, space infrared (5-210yum) tele¬ 
scope (target launch in 2022). With its 3.2-meter telescope cryogenically cooled to 6 
Kelvin, SPICA is 100 ti mes more sensitive than its precursors. Its Mid-infrared Camera 
and Spectrograph fMCS: lKataza et a02012l:IWada et alll2012h has a large field of view 
of 5x5 arcmin, and is sensitive in 5-38//m. In Figure |7J we compare SPICA/MCS’s 
survey speed (in 5x5 arcmin^ area) with that of James Webb Space Telescope (JWST). 
The figure shows that SPICA’s survey speed is comparable to JWST at < 15yum, and su¬ 
perior at > 18/rm. Note however at > 20yum, SPICA reaches galaxy confusion limit (the 
blue filled circles) in less than one hour. Taking advantage of the wide-held of view, 
and superb sensitivity of the MCS, we propose a wide area confusion-limited imaging 
survey using all the broad hlters of MCS. 
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2.1. Survey design 

2.1.1. Deep survey 


It has been known that infrared (IR) properties of galaxies depend on galaxy environ¬ 
ment. For example, iGotol ( 2005 ) showed that galaxy density distribution depends on IR 
luminosity, in a se nse that more IR lu minous galaxies are in less dense environment in 
the local Universe. iGoto et al.l(l2010a) showed that t he sha pe of the restframe 8/im LFs 
depends on the local galaxy density. Kovama et all (l2008h found the excess of 15-/rm- 
detected galaxies in the medium-density environments. These examples show that we 
need to investigate various environments from dense cluster cores to rarefied field, fo 
fully undersfand infrared properfies of galaxies and fheir evolufion. As a minimum of 
such sfudy field, we propose one deg^ survey area. Such a field covers 30 by 30 Mpc af 
z=3, including one massive clusfer in fhe survey area, allowing us fo invesfigafe from 
dense clusfer core fo fhe rarefied field continuously. 

This survey is also where fhe MGS is fechnically advanfageous; fhe MCS’s field 
of view (5 by 5’) is much larger fhan fhaf of fhe JSWT/MIRI (1.7 by 1.3’). Therefore, 
in ferms of fhe survey speed of 5x5’ area or larger, fhe MCS is comparably sensifive fo 
fhe JWST/MIRI even af <20//m (Fig. |7]). Af longer wavelengfh of >20//m, fhe MCS is 
more sensitive fhan fhe JWST/MIRI. 

We need fhe depfh of ~1 /rJy af 10 jum fo invesfigafe ULIRGs up fo z=3 . In order 
fo fully sample mid-IR SED, we aim fo image fhe area in all 8 broad fillers of WFC-S 
channel and all 7 broad fillers of WFC-L channel. Since in fhe WFC-L wavelengfh 
range, we reach fhe galaxy confusion limif in a shorf exposure lime (~15 min), and 
because fhe WFC-S and WFC-L channel can observe simullaneously, fhe imaging in 
fhe 7 fillers of fhe WFC-L channel can be oblained during fhe longer exposure of fhe 
WLC-S channel. In lolal, we need 9.67 hours per LoV (nel exposure lime of 1 hour per 
filler). We cover 0.9 deg^ by 12x12 tiling. This will fake 1392 hours, or 60 days. This 
is 6% of 2.5 years of SPICA mission lifetime. 

The survey needs fo be coordinaled fo be in fhe same region as fhe far-infrared 
SALARI survey. The SALARI survey needs ~1000 hours fo observe 1 deg^ region, fo 
fhe depfh of z=3 ULIRGs. Therefore, in lolal, we need ~2600 hours, or 3.6 monlhs of 
SPICA lime fo observe fhe field. This is a large amounl (12%) of fhe SPICA mission 
lifetime, and Iherefore, fhe field needs fo be of a good visibilily. Ancillary dala from 
olher wavelenglhs are desirable. We propose fhe north ecliplic pole or soulh eclipfic 
pole regions as candidales. IR cirrus confusion needs fo be smaller fhan galaxy confu¬ 
sion. We need fo avoid regions wilh >3 MJy/sr. 

The near-infrared LPC camera (5x5’ field of view, 5 broad fillers in 0.7-5 gm) can 
observe while fhe MCS is faking image (one of fhe HRS or MRS cameras need fo be 
off). Therefore, in Ihis region, near-IR image in 5 broad fillers of LPC can be oblained 
simullaneously. The depfh in one hour LPC exposure is expecled fo defecl ULIRGs af 
z=3. These near-IR pholomelry will be useful fo compute pholomelric redshifl of MCS 
sources combined wilh fhe ground-based optical imaging. 

The WLC-L will expose for one hour in each filler, while if reaches confusion 
limit in 15 min. Therefore, >4 times redundant data will be obtained in each filter of 
WLC-L. The survey needs to be carefully planned to sample variability of the sources 
in various time scales. 
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The survey volume at z=3 is 1.5x10^ Mpc^ (with Az of 0.2). The number density 
of ULIRGs at z=3 is 10“^ Mpc“^ dex“^ Therefore, we expect to detect 150 ULIRGs 
at z=3 in a slice of Az of 0.2. 

At z=l, the survey volume is 4.0x10^ Mpc^ (with Az of 0.2). The number density 
of ULIRGs measured by AKARI at z=l is 10“^ Mpc~^ dex“^ (Figl5l lGoto et al.ll2010bh . 
Therefore, we expect to detect 40 ULIRGs at z=I in a slice of Az of 0.2. 

2.1.2. Wide survey 

Complementary to the deep survey, the wide survey aims to detect rare objects that 
could not be found in numbers in the deep survey. For example, the number density of 
hyper luminous infrared galaxies at z=3 is 10“^ Mpc“^ dex“^. In the volume of deep 
survey (L5xl0^Mpc^, Az of 0.2 at z=3), there will be only 15 hyper LIRGs, and thus, 
not enough for statistical study. In the wide survey, we aim to survey 10 deg^ to the 
limit of 3-30 pJy. This will give us a 10 times larger volume, and thus, we will create a 
large enough sample for statistical study with 150 hyper-LIRGs. 

To reach 3 pJy, the total time requred is 0.975 hrs per FoV. 38 by 38 tiling can 
cover 9.05 deg^. Total telescope time required is then 1407 hours, or 6% of 2.5 years 
of mission lifetime. 

It is ideal that central part of the wide field to be covered by the deep survey, so 
that we understand the location of the deep field in terms of larger scale. However, 
both surveys require significant fraction (6% each) of mission lifetime, and both re¬ 
quire SAFARI counterpart surveys and various follow-up spectroscopic observation. 
Considering these factors, the deep and wide surveys are likely to be separated one at 
NEP and another near SEP. 
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